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O 

■ Masses of the p-wave excited heavy baryons have been calculated to the Aqcd /mg order using 

QCD sum rule method within the framework of heavy quark effective theory. Numerical results for 
kinetic energy Ai and chromo-magnetic interaction A2 are presented. The splitting between spin 
1/2 and 3/2 doublet derived from our calculation is given, for which the agreement with the current 
experiment is desirable. 
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INTRODUCTION 
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In the past decade continuous progress has been made in the investigation of excited heavy baryons. 
The lowest lying orbitally excited charmed states A c (2593) and A c (2625) have been observed by several 
collaborations lj, the excited states of E c have also been reported recently^, and more data are expected 
in the near future. From the theoretical prospect of view those data need to be studied comprehensively. 
Furthermore, with the collection of more experimental data for exited heavy baryon states it is useful to 
make some theoretical predictions on their spectroscopies. 



Heavy baryons containing a single heavy quark can be described exactly by heavy quark effective theory 
(HQET) 00, 0] in the heavy quark limit. This fact should be contributed to the spin-flavor symmetr 
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of system comprised of infinitely heavy quarks. HQET has been applied successfully to learn about the 
properties of heavy mesons and baryons, including the spectroscopy and weak decays. The mass formula 
for a spin symmetry doublet of heavy baryons up to order 1 /toq corrections can be written as 

M = m Q + A - (Ai + d M A 2 ) , (1) 

where parameter A is the effective mass of the light degrees of freedom in the wiq — > 00 limit, Ai and A2 
are related to the heavy quark kinetic energy and the chromomagnetic energy of HQET, respectively 

Ai = (B(v) I h v (iD x ) 2 h v I B(v)) , 
d M A 2 - (B(v) I h v h v I B{v)) . (2) 

The constant du characterizes the spin-orbit interaction of the heavy quark and the gluon field, it is zero 
for singlet and 1 , — | for spin 1/2, spin 3/2 doublet, respectively. Thus the splitting of the spin 1/2 and 
3/2 doublets is 

M%- Q ~ K Q = f A a > (3) 

where Bq, Bq denote spin 1/2 and 3/2 doublet, respectively. 

The heavy baryon mass parameters Ai and A2 play a most significant role in our understanding of the 
spectroscopy ;£>] and inclusive decay rates • They must be estimated in some nonperturbative approaches 
due to the asymptotic freedom property of QCD. A viable approach is the QCD sum rules 0] formulated 
in the framework of HQET 0|. This method allows us to relate hadronic observable to QCD parameters 
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via the operator product expansion (OPE) of the correlator. In the case of ground state heavy baryon, 
predictions on the mass spectroscopy have been made to leading and next-to-leading order in a s |lOl 
and to order l/mq 0,0] using QCD sum rule method. As to the excited baryon mass spectroscopy, 
only results to leading order in I/tuq expansion have been obtained from QCD sum rule HQ- In Q 
we have calculated the baryonic parameters Ai and A2 for the ground state Aq and Eq baryons using 
QCD sum rules in the HQET. Employing the baryonic currents from we now derive these parameters 
for excited A- and E-type baryons following the same procedure. 

The remainder of this paper is organized as follows. In Sec. Ill Al we introduce the interpolating currents 
for excited state heavy baryons and briefly present the two-point sum rules. The direct Laplace sum rules 
analysis for the matrix elements is presented in Sec. Ill Bl The Sec. Illll is devoted to numerical results and 
our conclusions. Some comments are also available in Sec. IIIII 



II. DERIVATION OF THE SUM RULES 
A. Heavy baryonic currents and two-point sum rules 

In this work we adopt those currents constructed from Bethe-Salpeter equation in Ref.[l4| 
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in which C is the charge conjugation matrix, r is the flavor matrix which is antisymmetric for A- type 
baryon and symmetric for E-type baryon, T's are some gamma matrices, and a, b, c denote the color 
indices. V can be chosen co-variantly as 

r' = 7^75 (5) 
for 'Egki, A Qkl , Aqki^qki doublets' spin 1/2 baryon, and 

r ' = r^ipi = -^(swi + itntitpx), (6) 

for Eq/ji, AqjcI, Aqkx, Eqki doublets' spin 3/2 baryon, in which gt PlPl and 7t A , 1 are perpendicular to 
heavy quark velocity v, defined as g tjiv = g pu - v p v v ^ tp: = 7^ - jv p . For the singlets A Qfe0 and T, QKQ 
the r' is simply unit matrix /. Notations used here to describe excited state heavy baryons are the same 
as those used in Ref. 0, Q|: k denotes Ik = 1 and Ik = whereas K denotes Ik — and Ik = 1, 
in which Ik is the orbital angular momentum describes the relative motion of the two light quarks and 
the orbital angular momentum Ik describes orbital motion of the center of mass of the two light quarks 
relative to the heavy quark. Q denotes heavy quark and the number in subscript is the total angular 
momentum of the light diquark system. As that for ground state baryons, the flavor configuration of A 
type baryon is antisymmetric and E type is symmetric. Depending on the number of derivatives or the 
form of r', interpolating currents can have forms different from those listed above and may also be used 
in applications [IH . 

In the following analysis we would use those currents to interpolate excited heavy baryon states. At 
the leading order of the 1 / uiq expansion they do not mix with each other even with the same quantum 
number, but to the next-to-leading order the mixing of interpolating currents will appear. In our subsequent 
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calculations we would only use those currents and did not consider the effect resulting from the mixing of 
interpolating currents for references [ToL 0, 0] have shown that the stability criterion for QCD sum rule 
applications excludes the existence of interpolating currents mixing or though there does exist the mixing 
the numerical result will not change drastically compared with the case without mixing. 

The baryonic coupling constant in HQET are also needed in our calculation, they are defined in form 
as follows 

(0\j\B(v))=Fu, (7) 

where | B(v)} denotes excited baryon state and u can be the ordinary spinor u or the Rarita-Schwinger 
spinor u a in the HQET corresponding to spin 1/2 or spin 3/2 doublet, respectively. Irrespective of an 
irrelevant constant factor in the leading order, the coupling constants for spin 1/2 and spin 3/2 doublet 
are equivalent for their identical spin-parity of the light degrees of freedom. 

In order to determine the effective mass of the excited baryons, we analyze the two-point correlator 
defined as 

d 4 xe lk -*(0 | T{j(x)j(0)} | 0) = l -±t T r[TT+]Il{u), (8) 

where k is the residual momentum and u> = v • k. For large negative value of u> II(a;) can be expressed in 
terms of perturbative and nonperturbative contributions. The nonperturbative effects can be accounted 
for by including quark and gluon condensates ordered by increasing dimension, which are the series of 
power corrections in the "small" \/uj variable. The Borel transformation in the variable u> can help to 
improve the convergence of these nonperturbative series. 

With those interpolating currents listed in Eq. Q it is straightforward to obtain the two-point sum 
rules: 

% k0 e" AA - /T = ^ M"c/T) + ^ (^G 2 ) + § e"# , (9a) 



Fl Qkl e- A — /T = ^- Sr(. c /T) - ±- 2 + ^ { a q ) 2 e-«* , (9c) 

Fl QK1 e^' T = 2 ™L 6r(. c /T) - £ ( ^ G2) + m£ {mf e -# f (9d) 

Fi QK1 e~^ T = ™£ S 7 (u c /T) - EL ( — G 2 ) + f {m? e"^ , (9e) 
y 7r 4 2tt z tt 2 

Fl QK0 = ^ ^/T) ~f,(^) + ^ {qq f . (9f ) 

V ku n l 8 7T 7T 4 

In calculations we adopted the gaussian ansatz for the nonlocal quark condensate to get the dimension 6 
condensate contribution. Dimension D > 6 condensates are not included. The functions 6 n (u) c /T) arise 
from the continuum subtraction and are defined in llT 



B. Sum rules for Ai and A2 

For the evaluation of the matrix elements Ai and A2 we consider the three-point correlators as follows 
i 2 ( d 4 x ( d 4 ye lk - x ^ k '-y(0 I T{j(x) h v (iD x ) 2 h v (0) J(y)} | 0} = I±^TY[tt+] T K (u,u/) , 



i 2 j d 4 x j d 4 ye tk - x - lk '-y(0 I T{j(x) h v G» v h v (0) j(y)} | 0) = d M ^Tr[TT+] T s (cj, w') , (10) 
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where the coefficients Tk {lo, lo') and Ts(lo, lo') are analytic functions in the "off-shell energies" lo — v-k and 
u) 1 = v ■ k! with discontinuities for positive values of these variables. Saturating the three-point functions 
with complete set of baryon states, one can isolate the part of interest, the contribution of the lowest-lying 
baryon states associated with the heavy-light currents, as one having poles in both the variables to and lo' 
at the value to = ui' = A. 

Confining us to take into account these leading contributions of perturbation and the operators with 
dimension D < 6 in OPE, the relevant diagrams in our theoretical calculations for the kinetic energy are 
shown in Fig. 1. The relevant diagrams for the chromo-magnetic interaction do not differ from those for the 
ground state baryons [l3| . so we do not show them here. Using dispersion relations Tk{lo, lo') and Ts(lo . lo' ) 
can be casted into the form of integrals of the double spectral densities. Following Refs. [Til UvL Eof . 
introduce new variables lo + — ^(lo + lo') and lj_ = lo — lo', perform the integral over w_, assume quark- 
hadron duality in w+, and employ Borel transformation B^, to suppress the continuum contributions 
and subtractions, we then obtained the desired sum rules. Considered the symmetry of the correlator it 
is natural to set the parameters r, r' to be the same and equal to 2T, where T is the Borel parameter of 
the two-point functions. We ended up with the set of sum rules 

A^^-A/t = ^3*.^ _ % { ^) + m * e - J, , (lla) 

#>~F*e-W = ^T-M-c/T) - ^G 2 > + <^> 2 e"^ , (lib) 

7T 7T 7T 7T 07T 



(lie) 



7T* 7T 2ir z IT 07T 

^e-^ = 2 ^^%(lo c /T) - ^G 2 > + {qq ? e-iS , (nd) 

_<«> F v^ . ^ M + f?<^> + M w . , (lle) 

7T 4 2 7T Z 7T lb 

7T* 4 7T Z 7T lb 



(llf) 



_ A ^Kx F2e -A/T = 2 S 9 (lo c /T) - ^(^G 2 ) + ^ fe) 2 e"^ , (llg) 

-X^F^ = ™^ 9 K/T) + |^(- G 2 > + *g { , q f e~£ , (llh) 

7T 4 8 7T Z 7T 64 

_ A A Qfc0jF2e -A/T = *^T^ ( /y) _ ±^<^G 2 ) + ^ ( g -,) 2 , (Hi) 

7T 4 8 7T Z 7T 64 

_ x ^ F 2 e -x,T = 26 32 nrl ° j 9 (^/r) - |^G 2 ) + (,-g) 2 e -# . (llj) 

7T 4 2 7T Z 7T 32 

The unitary normalization of flavor matrix Tt[tt + ] = 1 has been applied to get those sum rules, as what 
has been done for the two-point sum rules. 



III. NUMERICAL RESULTS AND CONCLUSIONS 

In the following analysis, the standard value for the condensates are adopted |8j. From two-point sum 
rules the effective mass can be obtained via a derivative of the Borel parameter as A = T 2 d\nE/dT, where 
E denotes the right hand side of the obtained two-point sum rule. Then comply to the standard procedure 
of sum rule analysis, we changed the continuum threshold lo c and Borel parameter T to find the optimal 
stability window and the numerical value of the effective mass will be determined within this window. For 
those sum rules obtained above, we found the typical value for the continuum threshold is lo c ~ 1.6 GeV 
and typical interval for the Borel parameter is AT ~ 0.4 GeV, which is a narrower one than the window 
for the ground state baryon. The only exception of this assertion is Aq/jq, for which we found a much 
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lower continuum threshold w c ~ 1.2 GeV and a narrower interval AT ~ 0.3 GeV. Also it is worth noting 
that for Agfco, the main contribution does not come from the perturbative part. The condensate ones play 
an important role in the determination of stability window. Indeed, if one keeps only gluon condensate 
contribution, there will be no stability window at all, just like the case in Ref. |l4) . But if assuming the 
perturbative dominance and omitting condensate contributions, we only ended up with a better stability 
and the numerical value is almost exactly the same. For the other sum rules for the effective mass the 
case is different. The dominant contribution to those sum rules comes from the perturbative part, and 
the dimension 4 and dimension 6 operators in the OPE only amount to 20% and 10% within the stability 
windows, respectively. The numerical results are presented in Fig. 3. For the aim of clarity we give our 
numerical average for the effective masses in Tabic [|] 

In order to get the numerical results for two 1/niQ order parameters, we divide our three-point sum 
rules by two-point functions to obtain Ai and A2 as functions of the continuum threshold oj c and the Borel 
parameter T . This procedure can eliminate the systematic uncertainties and cancel the parameter A. From 
the experiences of QCD sum rule applications in the field of heavy quark physics it is well known that 
three-point sum rule receives heavier contamination from the continuum modes than two-point one, and 
the stability is not as good as that for two-point sum rule El El 0. Our results for Ai sum rules of 
those excited heavy baryons are shown in Fig. 3 and in Fig. 4 results for A2 sum rules are presented. The 
typical value of continuum threshold is lj c ~ 1.8 GeV, io c ~ 2.4 GeV for Ai and A2 sum rules, respectively. 
And the typical interval of the stability window is also AT ~ 0.3 GeV for both. The Ai sum rule for the 
Y<qkq baryon is an exception of above statement, for which the typical value of continuum threshold is 
lo c ~ 1.4 GeV and the interval is AT ~ 0.2 GeV. As for the convergence of the OPE within the stability 
windows for those sum rules we would like to state some facts. For the Ai sum rules, except the case 
of AQfci, in which the dimension 4 and 6 operators give rise to contributions almost equal to that of the 
perturbative part, all other sum rules behave well, in which the contributions of the dimension 4 and 6 
operators amount to 50% and 20% of that of the perturbative part. On the other hand, for the A2 sum 
rules, the dimension 4 operator still contributes to almost 50% compared to that of the perturbative part, 
but the dimension 6 operator gives only a negligible contribution, typically less than 10%. The numerical 
results are listed in Table [I] too. 
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TABLE I: Effective mass in GeV, kinetic energy and chromo-magnetic interaction energy in GeV 2 for excited 
heavy baryons. Errors quoted are due to the variation of the Borel parameter T and continuum threshold ui c . 

With those values and heavy quark masses given in E3]> m c = 1-41 GeV and m& = 4.77 GeV, masses 
of excited heavy baryons to order 1/ttjq can be obtained immediately. We give those masses in Table ITT1 
Our results are comparable to the prediction on the excited heavy baryon masses obtained by using quark 
potential model |22^ . 
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6.199 
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TABLE II: Masses in GeV for excited heavy baryons. 



The splitting between spin 1/2 and spin 3/2 doublet can be obtained by multiplying A2 by a factor 
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3/2, cf. Eq. ©. It is 0.13 ± 0.02 GeV 2 , 0.16 ± 0.05 GeV 2 , 0.20 ± 0.02 GeV 2 and 0.19 ± 0.04 GeV 2 for 
Aqki, Aq^i, T^qki and Sgfei doublets, respectively. The approximately equal value for K.qki and Aq^i, 
Sqki and Sq^i may be interpreted as the signal which implies that current mixing effect cannot be large. 
If taking the middle value as theoretical predictions for the physical state, then the splitting for excited 
baryon states are 

a q 2 i " a qi = 0-15 ± 0.03 GeV 2 , 

Sq 2 - T? Q1 = 0.20 ± 0.03 GeV 2 . (12) 

Based on the current experimental dataji^], the splitting of excited A c doublet is 0.17 GeV 2 , which is in 
agreement with our theoretical result. When it is scaled up to the bottom quark mass scale there will be 
a factor ~ 0.8 approximately due to the renormalization group improvement. 

To conclude, we have calculated the 1/toq order corrections to the excited heavy baryon masses from 
QCD sum rules within the framework of the HQET. From thus obtained spectrum for the c quark case, we 
found that Aq^o and Eq/ci baryons lie ~ 600 MeV above the ground state baryon A c , while Aq^i, Aq^i, 
^QKi he ~ 800 MeV above A c and for T,qki baryon this value is ~ 900 MeV, typically with an error 
~ 300 MeV. When it comes to the b quark case, the result is that Aqw and Sq^i baryons lie ~ 300 MeV 
above the ground state baryon A&, while Aq£i, Aqki, ^qki and ^qki he ~ 550 MeV above Aj, for which 
the typical error is ~ 200 MeV. For the c quark case, 1/mg order corrections are ~ 400 MeV, which is 
not a small one due to the large value of the kinetic energy. When it comes to the b quark case those 
corrections will be suppressed by the still larger b quark mass. Our theoretical predictions for the doublet 
splitting are in agreement with the current experimental data. 
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Figure Captions 

Fig. 1. Non- vanishing diagrams for the kinetic energy. The kinetic energy 
operator is denoted by a white square, the interpolating baryon currents by 
black circles. Heavy-quark propagators are drawn as double lines. 

Fig. 2. Sum rules of the effective mass A for: (a) Ag^o, (b) Ag^i, (c) Eg/a, 
(d) Ag^i, (e) Sgifi and (f) Sg^o baryons. The different choices of the contin- 
uum threshold oj c corresponding to different curves are designated in individual 
figures respectively. Curves are plotted against the Borel parameter T. 

Fig. 3. Sum rules of the kinetic energy for: (a) K Qk0 , (b) Ag fe i, (c) Egfci, (d) 
Agifi, (e) Yjqki and (f) T,qko baryons. Others are the same as those in Fig. 
2. 

Fig. 4. Sum rules of the chromo-magnetic interaction for: (a) Ag^i, (b) Egfci, 
(c) Kqki and (d) T,qki baryons. Others are the same as those in Fig. 2. 
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